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Abstract 

Purpose The purpose of this paper is to enhance the math¬ 
ematical and physical understanding of practitioners of un¬ 
certainty analysis of life cycle inventory (LCI), on the 
application of possibility theory. The main questions dealt 
with are (1) clear definition of the terms—“necessity- 
possibility,” “probability,” “belief-plausibility,” and of their 
mutual relationships; (2) what justifies the substitution of 
classical probability for possibility; (3) mutual comparison 
of, and transformations in both senses between probability 
and possibility uncertainty measures; (4) how to construct 
meaningful input possibility measures from available prob¬ 
abilistic/statistic information; and (5) comparative analysis 
of the solutions of the problem of data uncertainty propaga¬ 
tion in LCI, afforded, respectively, by probabilistic Monte 
Carlo simulation and possibilistic fuzzy interval arithmetic. 
Methods The questions above are addressed from the rigor¬ 
ous mathematical formulations of the theories of probability 
and statistics, of possibility, and of random sets and belief/ 
plausibility functions, although directed to LCI uncertainty 
analysis practitioners. On this respect, the paper allows two 
different levels of reading: a basic level (main text) and a 
deeper level (Electronic supplementary material). Particular 
tools used are (a) various transformations between possibil¬ 
ity and probability distributions, in both senses, for the 
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continuous case, proposed by Dubois et al. (e.g., Reliable 
Comp 10:273-297, 2004); (b) Monte Carlo simulation for 
either independent or dependent input random variables; (c) 
fuzzy interval arithmetic; and (d) Heijungs and Suh’s (2002) 
matrix formulation of LCI problems. 

Results and discussion The links among uncertainty meas¬ 
ures, uncertain variables, and uncertainty analysis are 
cleared up. It is recalled how a probability measure can be 
constructed and attached to an input variable, and its prob¬ 
ability distribution and unknown “correct value” be related, 
in a physically meaningful way. It is justified that, usually, a 
dual necessity-possibility measure has much less uncertain¬ 
ty information than a comparable probability measure. Al¬ 
though the specialists are not unanimous, it is opined that 
the theoretical framework developed by Dubois et al. (e.g., 
Reliable Comp 10:273-297, 2004) is the most convenient 
one to use in uncertainty analysis, to compare and mutually 
transform probability and possibility data. This is exempli¬ 
fied in (a) the transformation of the very common triangular 
possibility and normal standard probability distributions; (b) 
the general construction of possibility measures from differ¬ 
ent probability data previously available; and, above all, (c) 
the comparison of the output information of possibilistic and 
probabilistic uncertainty analyses of an LCI problem pro¬ 
posed by Tan (Int J Life Cycle Assess 13:585-592, 2008). 
The general problem of data uncertainty propagation 
through deterministic models (e.g., of LCI) is tackled with 
(1) classical probabilistic Monte Carlo simulation (for either 
independent or dependent input random variables); (2) pos¬ 
sibilistic fuzzy interval arithmetic; and (3) hybrid methods 
(only mentioned). 

Conclusions (1) The practical conditions in which an anal¬ 
ysis of uncertainty should switch from a probability to a 
possibility basis are still ill-defined, but that seems to be the 
case when the input information is based on states of large 
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ignorance. (2) A dual necessity-possibility uncertainty mea¬ 
sure can be viewed both as an imprecise probability measure 
that substitutes a definite probability for an interval, and as a 
belief-plausibility measure. (3) A possibility distribution 
can be mathematically and physically interpreted as a ran¬ 
dom set of nested prediction/confidence intervals for the 
“correct value” of the variable, with confidence levels rang¬ 
ing from 0 to 1. (4) There exist mathematically and physi¬ 
cally sound rules to compare/transform probability and 
possibility uncertainty information under different applica¬ 
ble paradigms (e.g., based on the reliability of the input 
information). (5) Sometimes, Geer and Klir’s (Int J Gen 
Syst 20:143-176, 1992) confidence index has a physically 
counterintuitive behavior in uncertainty analysis. (6) The 
probabilistic Monte Carlo simulation can be used also for 
dependent random input variables, only requiring more ex¬ 
igent input information (conditional probability distribu¬ 
tions) than in the case of independency. (7) The 
possibilistic fuzzy interval arithmetic uncertainty analysis, 
although computationally cheap, generates output informa¬ 
tion quite poor—a good point estimate but a set of (roughly) 
confidence intervals with very large amplitudes for the 
“correct value” of each output variable. 

Recommendations (1) In probability uncertainty analysis, 
pay attention to the relation between the “correct value” of 
a random variable and the parameters of its probability 
distribution (e.g., mean or mode). (2) Do not precipitate in 
changing from probabilistic to possibilistic uncertainty anal¬ 
ysis: it may be theoretically unjustified and much output 
uncertainty information is lost. (3) Respect the well- 
established applicable rules in going from probability to 
possibility uncertainty information, or vice versa. (4) Be 
attentive to possible counterintuitive physical meaning of 
Geer and Klir’s (Int J Gen Syst 20:143-176, 1992) confi¬ 
dence index in possibilistic uncertainty analysis. 

Keywords Belief-plausibility measures • Fuzzy • Interval 
arithmetic • Monte Carlo simulation • Necessity-possibility 
measures • Probability-possibility transformations 

1 Introduction 

The need to or the interest in performing uncertainty analy¬ 
sis of life cycle inventory (LCI) is well known and does not 
require justification at this time (e.g., Lloyd and Ries 2007; 
Reap et al. 2008a, b). 

The most informative tools of data uncertainty analysis 
supply interval estimates for the “correct values” of the 
output variables of interest, with quantitative measures of 
uncertainty attached. Classically, the “uncertainty measure” 
chosen is a “probability,” by building up the uncertainty 
analysis upon the mathematical theory of probability and 


statistics (e.g., Peters 2007). However, recently, an alterna¬ 
tive “possibility measure” is being proposed by some 
authors (e.g., Geldermann et al. 2000; Tan et al. 2002; 
Benetto et al. 2006; Tan 2008), which is defined in the 
theory of possibility, itself a relatively modern mathematical 
theory (Zadeh 1978; Dubois and Prade 1988), with clear 
though subtle and easily misleading relations with other 
modem mathematical theories of fuzzy sets (Zadeh 1965), 
random sets and belief/plausibility functions (Shafer 1976), 
and, of course, the old theory of probability and statistics. 
This paper intends to enhance the mathematical understand¬ 
ing and contextual physical interpretation of these theories 
in uncertainty analysis of LCI. It is directed both to soft 
(main text) and deep (see Electronic supplementary materi¬ 
al) practitioners. The authors feel particularly indebted with 
Dubois et al., whom publications in this field cover a prac¬ 
tically continuous period of 30 years up to the present. In 
this respect, Dubois et al. (2000) deserves to be highlighted. 

With the above purpose, the paper commences (Section 2) 
with the presentation of the general backstage LCI problem 
and, within it, of a specific backstage case. In Section 3, the 
basic relations among uncertainty measures, uncertain vari¬ 
ables, and uncertainty analysis are set up, and an overview 
of subsequent Sections 4 to 7, through the formulation of the 
two main problems of uncertainty analysis dealt with, is 
given. In Section 4, it is recalled how a probability measure 
can be constmcted and related with the distribution of a 
random variable, and what is the relationship between the 
latter and the “correct value” of the underlying physical 
variable. In Section 5, dual possibility-necessity measures 
are introduced, and it is detailed both how they can be 
constmcted for input variables (Section 5.1) and how they 
should be interpreted for output variables in uncertainty 
analysis (Section 5.2). In Section 6, the key question about 
which uncertainty measure, probability or possibility, can or 
should be attached to a given input physical variable is 
answered. In Section 7, the following issues of the main 
problem of propagation of uncertainty from input to output 
variables are dealt with (1) the possibility approach of fuzzy 
interval arithmetic (Sections 7.1 and 7.4); (2) the probability 
approach of the simulation of Monte Carlo (Sections 7.2 and 
7.4); (3) the comparison of both former approaches in the 
uncertainty analysis of the backstage case (Section 7.3); and 
(4) emerging mixed probability/possibility approaches (Sec¬ 
tion 7.5). In Section 8, the main conclusions are highlighted. 

2 Backstage LCI problem 

Let us put in the backstage the general LCI problem pertain¬ 
ing to an “economic system” (e.g., a plant unit) encompass¬ 
ing n “economic processes” that consume/produce n input/ 
output “economic commodities” (e.g., secondary materials, 
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products) while extracting/emitting from/to the ambient m 
“environmental commodities” (e.g., raw materials, pollu¬ 
tants), with the following matrix formulation (Heijungs 
and Suh 2002): 

g = B A 1 f (1) 

In the matrix Eq. l,A(n x n ) is the “technology matrix” of 
the system, with general element ay= (negative/positive 
amount of input/output economic commodity i consumed/ 
produced per unit “scaling level” or time of operation of the 
economic process J), here supposed to be invertible, which 
implies that the, so-called allocation problem does not arise 
or has already been solved out for the system; B(m x n) is the 
“intervention matrix” of the system, with general element 
b k j= (negative/positive amount of input/output environmen¬ 
tal commodity k extracted/sent from/to the environment per 
unit scaling level or time of operation of the economic 
process j)\ f{n x 1) is the “functional unit column-matrix” 
specified in the LCI problem, with general element f t = 
(negative/positive amount of input/output economic com¬ 
modity i that the system is demanded to consume/produce); 
and the output g(m x 1) is the “inventory or environmental 
interventions column-matrix,” with general element g k = (to¬ 
tal negative/positive amount of input/output environmental 
commodity k extracted/sent from/to the environment by the 
system). 

Tan’s (2008) case I will serve as a specific example 
(Fig. 1), heretofore named “backstage case.” In this case, 
the system encompasses processes I and II (j= 1 and 2). Both 
processes emit the pollutant D (environmental commodity 
with k= 1) to the environment, but while process I consumes 
the product B (economic commodity with i= 2) and produ¬ 
ces the product A (economic commodity with / = 1), process 
II does the inverse, which allows the system to operate in an 
economically efficient recycling mode. The technology A , 
intervention B , and functional / input matrices of the sys¬ 
tem/problem are given below (Eqs. 2a-c), with all amounts 
expressed in mass units of kilogram. For each one of the 
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Fig. 1 Scheme of the so-called, in the paper, backstage case: Tan’s 
(2008) case I system operating in the steady-state condition demanded. 
(Notice that for this condition to be physically possible, some inputs 
must be neglected, which does not matter here.) 


“uncertain input variables” (a n , a 22 , b n , Z? 12 ), best estimates 
of the lower bound, point value, and upper bound of the 
respective (unknown) “correct value” are given. The single 
scalar output variable is gi= (total amount of pollutant D 
emitted by the system for this to satisfy the required 
demand). 
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100 " 

0 


[(0.05,0.06,0.07) (0.9,1,1.1)], 


(2a-c) 


3 Uncertainty measures and uncertain variables 
in uncertainty analysis 

Any “uncertain variable” X with domain C/cR (R being the 
mathematical set of real numbers) can be associated with an 
actual or virtual “non-deterministic experiment” E with set 
of results (7, i.e., a “repeatable” experiment in which, al¬ 
though all relevant “deterministic” factors are directly or 
indirectly controlled, the value x&U taken by X in any 
individual realization of E (in uncertainty analysis, usually, 
x is the outcome of a measurement procedure specified in E) 
cannot be predicted deterministically in a reliable way (e.g., 
Andre 2008). In any realization of E, if its result is xgA then 
“event” A is said to have occurred, and M(v4)e[0, 1] is the 
general designation of a quantitative measure of the experi¬ 
menter’s (un)certainty that A will actually occur. For in¬ 
stance, for the variable X=a n of the backstage case 
(Section 2) and the event A = [0.5, 1], whatever measure M 
is adopted, one would expect that M(v4)=l. 

The first problem the uncertainty analysis tools dealt with 
in this paper must solve is to choose an appropriate uncer¬ 
tainty measure M connected with each uncertain variable X 
and respective non-deterministic experiment E with set of 
results U. Two main alternative solutions to this problem are 
considered here: either a probability measure P (Section 4) 
or a dual necessity-possibility measure N-II (Section 5), 
both defined in U. The choice between both is discussed 
in Section 6. 

The second problem the uncertainty analysis must solve 
is the propagation of uncertainty from input to output vari¬ 
ables, i.e., “Given: n uncertain input variables (X u X 2 , ..., 
X n ) with known uncertainty distributions, and an output 
variable Y related with the former input variables through 
a known deterministic function Y=h(X\, X 2 , ..., X n ); obtain: 
the uncertainty distribution of Z” Consistently with the first 
problem, some or all the input variables can be probabilistic 
or random variables and, the remaining ones, possibilistic or 
fuzzy variables. Besides, different probabilistic/possibilistic 
relations may exist among the input variables, ranging from 
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simple independence to different types and degrees of cor¬ 
relation. This problem is dealt with in Section 7. 

4 Probability uncertainty measures and random 
variables 

As a sound foundation of the mathematical theory of prob¬ 
ability and statistics, Kolmogorov (1950) defined a “proba¬ 
bility measure” P in a universe set U, , by associating to 
(almost) any subset A of U, a number P (A), satisfying the 
following three axioms: (Al) P(^4)>0, (A2) P(£/)=l, and 
(A3) (additivity axiom for two sets) if AC\B=0 , then P 
04u£)=P(yt)+P(5). 

To attach a probability measure P (in U) to a given 
uncertain physical variable X with domain UQ R, in a phys¬ 
ically meaningful way, it is required, in the first place, that X 
be a “random variable,” i.e., that the non-deterministic ex¬ 
periment E behind X falls into the restrict class of “random 
experiments.” 

In the second place, the measure of probability P linked 
to the random experiment E with set of results U can be 
constructed from scratch through mainly three methods 
(Andre 2008) —method 1 (Empirical Method): to perform 
E as many times as possible («—> 00 ) and to rely upon the 
asymptotic identity: P(A ) = lim f n (A)/n , where f n (A) = 

n—> 00 

(number of occurrences of A in the first n realizations of 
E); method 2 (Theoretical Method): if U is finite and the 
physical “hypothesis of symmetry” according to which “all 
results r in U are equally likely” is applicable, then P(^4)=n 
( A)/n(U ), where n(-)= “number of elements of;” and method 
3 (Subjective Method): P(v4)= (subjective estimate of an 
“expert”). It is important to understand that each one of these 
methods has its own non-mathematical conditions of correct 
application. For instance, in method 1, problems with the 
control of experiment E can be detected through the apparent 
divergence of the succession p n =f n (A)ln , and, in method 3, the 
choice of the “expert” can be very relevant. Besides, if all 
methods fail to give any or mutually consistent results, then, 
the possibility that the experiment E is intrinsically “non- 
random” should be considered (cf. Section 6). 

In the third place, the relation between the probability 
measure P attached to (E, U) and the probability function p\ 
t/—>[0,1] of the random variable X, if X is discrete, or the 
probability density function / of X, f.U—>[ 0,+oo[, if X is 
continuous, is, respectively: 

m = 5>w or / /(x) • dx (3a, b) 

xeA J. 

xeA 

Finally, in the fourth place, in uncertainty analysis, denot- 
ing by x the “correct value” of X , if its measurement 
procedure is “non-biased” (in average), then x =/x, p being 


the probabilistic mean of X. Furthermore, if the probability 
distribution of X is symmetrical and unimodal then it is also 
x =/i=x 0 , x 0 being the “mode” of X. Both the former con¬ 
ditions apply to many physical random variables. However, 
non-symmetrical lognormal probability distributions are as¬ 
sumed for the variables of all unit processes in ecoinvent 
data vl.l (Frischknecht et al. 2005). Moreover, when X is 
measured indirectly through other basic random variables 
(X u X 2 , ..., X m ) with which it bears the deterministic rela¬ 
tionship X=h(X u X 2 , X m ), if h is non-linear and the 
standard deviations of the basic variables are not all small, 
the distribution of X will be biased and asymmetrical, even 
if the distributions of the basic variables are non-biased and 
symmetrical, i.e., x =h(x 1, x 2 , x m )^/x, x 0 even if 
x* i =ti i =x 0i (/=1, 2 , ..., m). 

5 Possibility uncertainty measures and possibilistic 
or fuzzy variables 

In the theory of possibility (Schackle 1961; Zadeh 1978; 
Dubois and Prade 1988), a “possibility measure” n is de¬ 
fined in the universe set U by associating to any subset A c U 
a number n (A) satisfying the following axioms (Dubois et 
al. 2000): (Al') n(0)=O; (A2') n(£/)=l (normal case); and 
(A3') (maxitivity axiom) for any A, BqU , n(Aui?)=max 
{11(A), n(i?)} (cf. the additivity axiom A3 in Section 4). 

If UqR, the possibility measure n can be associated with 
a “possibilistic or fuzzy variable” X characterized by the 
“possibility distribution” 7 r :R—>[0,1], such that: 

n(A) = sup{;r(x),x G A }, (4) 

with 7r(x)=l for one (and only one, the mode x 0 , if 7 r is 
unimodal, as is most usual) or more values xeR. The possi¬ 
bility measure n has also associated to it the dual “necessity 
measure” N defined, in this way: 


N(A) = 1 - n(A) 

(5a) 

It is simple to show that (Dubois 2006): 


N{A) = inf {1 — n(x ), x € A } and 

(5b,c) 

N(A flS) = min{N(yl), N(5)} 



The identity Eq. 5c is called the minitivity axiom for the 
necessity measure (cf. the maxitivity axiom A3' above). 
Figure 2 illustrates the definition of the dual measure H-N 
for a generic unimodal possibility distribution 7 r. 

5.1 Construction of possibility measures for input variables 

Due to the relative youth of possibility theory (Zadeh 1978), 
the empirical procedures and physical hypotheses and ideas 
available to directly construct the possibility measure n 
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Fig. 2 Definition of the dual possibility-necessity measure II-N associated 
with the possibilistic or fuzzy variable X with unimodal possibility distri¬ 
bution 7 r with mode x 0 , for a generic event A. The event A a is called the “a- 
cut of 7r”. Per definition, irrespective of a , II (A a )= 1 and N(4 rv )= 1 -a 


attached to a physical input variable X are not so well 
understood as the ones behind the methods 1-3 that are 
used to construct probability measures (cf. Section 4). They 
are not reviewed here but the interested reader can consult, 
e.g., Hersh and Caramazza (1976), Hisdal (1991), and Tong 
et al. (2004). Instead, three useful indirect methods of con¬ 
struction of physically meaningful possibility measures are 
proposed below. In a more (methods I—II) or less (method 
III) extent, all methods start from a given probability mea¬ 
sure P (here X is still envisaged as a random variable) and 
afford, through an appropriate transformation n=T(P), the 
“nearest” possibility measure II (X is now seen as a possi¬ 
bilistic variable), in a certain sense. A basic hypothesis 
shared by the three methods, with the exception of the final 
variant of method III (triangular possibility distributions), 
consists in supposing that, for unimodal probability distri¬ 
butions P, the mode x 0 is the best point estimate x* of the 
“correct value” x of the physical variable behind X. In fact, 
if P is asymmetric then, most commonly, x =/i^x 0 (cf. 
Section 4, in fine). However, in these cases, it seems that a 
possible way around could be to translate the probability 
distribution P of X in order to enforce that xo = x*. The 
underlying theory is covered in Section A2.2 of the Elec¬ 
tronic supplementary material for the interested reader, 
where the rationale behind some symbology is also given. 

Method I A sound probability measure P is available for X , 
typically obtained with the method 1 (empirical) of Sec¬ 
tion 4. Here, we consider only the case of a continuous 
random variable X with unimodal probability density func¬ 
tion / with mode x 0 , strictly increasing for x<x 0 and strictly 
decreasing for x>x 0 . In this case, the possibility distribution 
7r of X should be obtained by applying to / the following 
transformation 7r=T 2 (f): 

n(x) = 7i(y) = F(x) + 1 - F(y) for x < x 0 ,y > x 0 and f{y) =/(x) 

(6a, b) 


In Eqs. 6a, b: F is the distribution function of variable X; 
and the interval [x, y\ can be interpreted as the/(x)-cut of f i.e., 
the analogue of the a-cut A a of i r (cf. Fig. 2). 

Method II A probability distribution P much less reliable 
than the one of method I, typically obtained with method 3 
(subjective) of Section 4, is available for X. Now, for a 
continuous random variable X with probability density func¬ 
tion f the transformation to apply is tt= 

7r(x) = / min{/(x),/(x / )} • dx' for any x E U (7) 

x'eU 

In Fig. 3, the possibility distribution tti resulting from the 
application of I) -1 to a normal standard random variable X 
can be compared with the possibility distribution tt 2 
obtained by method I for the same case. 


Method III The uncertainty information available on X is a 
set of reliable prediction intervals z*, (a'=a\ 9 ..., a' p ) of 
estimation (either empirical or subjective) of the “correct 
value” x of X , for some confidence levels a'<=[0,1]. In 
particular, for a-0 the prediction interval collapses into a 

simple point estimate: zj = 5*. If the prediction intervals (for 
cf>0) are classical statistics confidence intervals obtained 
from an empirical sample of X, for unimodal probability 
distributions P of X (in this case, unknown in detail), an 
estimator for the mode should be used and the sample must 
also be one-dimensional (i.e., have a single value). In any 
case (postulating also that 7r is unimodal with mode x 0 , 
increasing for x<x 0 and decreasing for x>x 0 ), the same 
transformation T 2 of method I can be applied to each pre¬ 


diction interval z*, = 
identities: 


Xn X i 


%M,a' 


, which leads to the simple 



— ^ (-^M, or') — 1 


a 


(8a, b) 



Fig. 3 Possibility distributions 7T2=T 2 (f) and iri = Ti~ l (f) obtained forX 
(envisaged as a possibilistic or fuzzy variable) from a previously 
available normal standard probability density function f of X (envis¬ 
aged as a random variable), respectively, by methods I and II. In both 

• • • • sfc 

cases, x o =0 is the best available point estimate of the “correct value” x 
of the physical variable behind X 
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In particular, for a'=0, Eqs. 8a, b further collapse into: 
/r(x*) = 1 (8c) 

The possibility distribution 7r(x) can be guessed in the rest 
of U with an appropriate best fitting procedure of the set of 
points {(x mt a',7r(x m ^)), (x M ,a',n(x M ,a')) for a'=a\, 
a'p). In the often case in which the only available informa¬ 
tion on X is v* and i\ = [x w j,xm,i] , it is a good option to 
take a triangular possibility distribution 7r for X, with mode 

xo = x* and support [x m , x M ]=i\ (e.g., the possibility distri¬ 
bution 7 r in Fig. 4), as justified in Section A3 of the Elec¬ 
tronic supplementary material. 

5.2 Physical interpretation of possibility measures of output 
variables 

Let us now suppose that the pertinent uncertainty propaga¬ 
tion problem formulated in general terms in Section 3 and 
dealt with in depth in Section 7 has already been solved, 
affording the possibility distribution n of the possibilistic 
output variable Y=h(X\, X 2 , ..., X n ). How should the uncer¬ 
tainty information inherent to tt be interpreted? To simplify 
without significant loss of practical generality, when not 
mentioned otherwise, 7r is supposed to be unimodal with 
mode y 0 . Besides, consistently with Section 5.1, it is pre¬ 
supposed that the mode x 0i of the possibility distribution 7iy 
of input variable X f (i= 1, 2, ..., n) is the best available point 
estimate of its “correct value” x*, i.e., x* = xo,. 

In the first place, the mode y 0 is the best available point 
estimate of the “correct value” y of the physical variable 
behind LJ namely, yo =y* =/z(x^,x^, ...,x^). 



Fig. 4 Geer and Klir’s (1992) “confidence index” C (A) (thick lines ) in 
events A=[xi,x 2 \c U, x\ and x 2 being, respectively, fixed and free 
parameters, for a fuzzy variable X with triangular possibility distribu¬ 
tion 7 r with support [0.5,1] and mode x o =0.6 (thin line). Positive/ 
negative values of C (A) are said to confirm/disconfirm event A but, in 
the case of interrupted lines, this interpretation of C(A) is counterintuitive 


In the second place, and this is the deepest applicable 
theoretical result, the a-cut A a of the possibility distribu¬ 
tion 7 r (see Fig. 2) can be roughly interpreted as a “con- 

fidence interval” /*, of estimation of y with confidence 
level a ,= l—a, as the latter is defined in classical statistics 
theory for a single-value sample. However, depending on 

7 r, the amplitude of /*, can be much larger than the 
amplitude of the optimum (i.e., narrowest) confidence 
interval with level a'. For instance, due to an optimal prop¬ 
erty of transformation T 2 , the a-cut of the possibility distribu¬ 
tion 7r 2 of A represented in Fig. 3 is the narrowest confidence 
interval for x with confidence level a ,= l~a, and, in this 
sense, 7r 2 has much more uncertainty information than tt \. The 
theory underlying this result is given in Sections A2.2 and A4 
of the Electronic supplementary material (see also: Yager 
1992, Smets 1990, Dubois et al. 2008). 

Finally, in the third place, the following “confidence 
index (in events)” C(^4)e[-l,l] of Geer and Klir (1992) is 
sometimes (e.g., Benetto et al. 2006) used to synthesize the 
uncertainty information contained in a dual possibility- 
necessity measure n-N: 

C(A) = 17(A) + N(A) - 1 = n(A) - n(A) (9a, b) 

Specifically, positive/negative values of C (A) are inter¬ 
preted as “degrees of confirmation/disconfirmation” of A. 
However, this interpretation of C(A) can sometimes be counter 
intuitive. For instance, in the case of a simple triangular 
possibility distribution 7 r with support U= [0.5,1] and mode 
x o =0.6 (see Fig. 4), for the family of events A= [0.53 ,x 2 ], C(A) 
keeps a constant value of 0.7 for increasing x 2 above 0.88, 
while, by physical sense (perhaps, induced by probabilistic 
thinking), one would expect that it continued to approach 1. In 
fact, the maximum degree of confirmation C(A)= 1 is only 
attained if, and only if, the events is impossible, with n (A)= 1 
and II(A)=0, while, for the former family of events A, no 
matter how large is x 2 , it is always n(v4)>7r(0.53)=0.3>0. 


6 The choice between probability and possibility 
uncertainty measures 

Let us commence by clearing up a basic question. Although 
Klir and Parvis (1992) (see Klir 2006 for a more recent 
presentation) say that (in a given universe set U) there exist 
measures of possibility-necessity n-N and of probability P 
that are “commensurate,” i.e., that have essentially the same 
uncertainty information, the authors follow Dubois et al., 
who defend convincingly that the former is always 
“weaker,” i.e., it has less precise uncertainty information, 
than a comparable probability measure (e.g., Dubois et al. 
1993, 2000; Dubois 2006). Their arguments are reviewed in 
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Section A1 of Electronic supplementary material (see also: 
Walley 1991, Giles 1982, De Cooman and Ayels 1999). 

The main question dealt with in this section is the following 
one: Given an input physical variable X associated with the 
non-deterministic experiment E, should the uncertainty inher- 
ent to the determination of its “correct value” x be quantified 
through a probability measure P or through a dual possibility- 
necessity II-N measure? In other words, should Xbe treated as 
a random variable or as a possibilistic/fuzzy variable? 

At the onset, two circumstances may dictate the conve¬ 
nience of preferring a possibility to a probability measure 
for X: (1) only method 3 of Section 4 (subjective estimates 
of experts) can be used to construct the probability measure 
and the available experts are largely ignorant on the out¬ 
comes of experiment E ; and (2) it is wished to reduce the 
computation time of resolution of the subsequent uncertain¬ 
ty propagation problem. Regarding case 1, notice that an 
expert totally ignorant about event A (in E) could consis¬ 
tently postulate that Ii{A)=Ii{A)=\ and N(A)=N(A )=0 (cf. 
Section 5) but if it postulates, for instance, that P(X)=0.8 
then it must consistently postulate that V{A)=\~V{A)=02 
(cf. axioms A1-3 of Section 4). Similarly, the so called 
“vacuous possibility measure,” according to which II (A)= 
1 for any AqU, is much less informative than the analogous 
“uniform probability measure” (Dubois 2006). In Section 
A1 (Electronic supplementary material), it is explained in 
more detail why possibility theory is a theoretical frame 
more adequate than probability and statistics theory to deal 
with uncertainty under states of large ignorance. Case 2 is 
discussed in Section 7.3. 

If the former circumstances simply do not occur or are 
not determinant, a random experiment E can be assumed 
and any of the methods 1-3 of Section 4 used to construct a 
probability measure P attached to X, or, what is equivalent, a 
probability distribution for X. In this process, two other 
circumstances, and none more the authors can think of, 
can still lead to the mandatory substitution of the probability 
model by a possibility model: (3) in the course of the 
application of method 1, E is suspected to be non-random, 
or (4) while applying method 3, independent estimates of 
the probabilities of a set of interconnected events (e.g., A , B , 
AP\B and AuB), supplied by one or more experts, lead to 
a gross violation of axioms A1-3 (Section 4). The pos¬ 
sibility of case 3 is already foreseen in Section 4. On 
case 4, Dubois (2006) mentions a study of Raufaste et al. 


(2003) that would suggest that subjective estimates of 
uncertainty measures tend to obey the maxitivity/mini- 
tivity axiom (A3') of possibility theory rather than the 
additivity axiom (A3) of probability theory. Taking the 
example given above within parenthesis, to postulate that 
P(X)=0.2, P(£)=0.4, P(XfTg)=0.1 but P(Xu5)=0.3 is 
inconsistent with axioms (A 1-3) because they imply that 
P(Xu5)=P(X)+P(£)-P(Xn£). 


7 Uncertainty propagation problem 


7.1 Possibilistic approach 


Let us take once again the backstage case (cf. scheme Fig. 1 
and numerical inputs Eqs. 2a-c). For each uncertain input 
variable X, following the last variant of method III of Sec¬ 
tion 5.1, Tan (2008) adopts a triangular possibility distribu¬ 
tion 7r with support [x m , Xm\ and mode x 0 , where ( x m , x 0 , Xm) 
are the best available estimates, respectively, of the lower 
bound, point value and upper bound of the “correct value” 
x of X(e.g., in Fig. 4 it is shown the possibility distribution 
7r postulated for X=an[x m =0.5, x o =0.6, xyv/= 1.0]). 

For reasons and in conditions (really, still not fully under¬ 
stood) that are explained in Section A5 of the Electronic 
supplementary material (see also: Zadeh 1975a,b, Gebhardt 
and Kruse 1994), the solution of the uncertainty propagation 
problem Y=h(Xi, X 2 , ..., X n ) in possibility theory can be cast 
into a simple form called after “fuzzy interval arithmetic,” 
used, namely, by Dubois et al. (2004a, b) and Tan (2008), 
which affords any a-cut A Y>a (<xe[0,l]) of the possibility dis¬ 
tribution 7r of the output variable k in the following form: 


A 


Y,a — 


y'aX'a 


with 


(10a) 


y a , y' a = inf, sup{/;(xi, G A x , a n ... f]x n G A x „, a } 

(10b, c) 

In Eq. 10b,c, Ax h a denotes the a-cut of the possibility 
distribution 7r z - of the input variable X t . Specifically for the 
general backstage LCI problem with matrix-Eq. 1 , in which 
Y=g k , the following expressions can be used to explicitly 
solve the min-max problem posed by Eq. 10b,c: 
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with : 


dgk 

day 


-Ckt ■ Sj and 


dgk 

dby 



(12a, b) 


In the right members of the multiple Eq. 11 : h stands for 
the scalar deterministic function linking the output g k to all 
the inputs (ay, by,/), consistent with the matrix-Eq. 1; and it 
is supposed that the a-cuts of the possibility distributions of 
the output variable g k and of the input variables ay and by 

are, respectively. Ag k ^ a gk,a ->&k,a ^ij,a i^ij,a 


and A bjJt a = [%,/,%(/'] • In Eq. 12a,b, which are easy to 
derive from Eq. 1 (Heijungs and Suh 2002), c ki and Sy are the 
general elements of matrices, respectively, C = B • A~ ! and 
s = A -1 • f, the latter also called “scaling or times column- 
matrix,” and S ki is Kronecker’s delta, i.e., S ki =0 or 1 whether 
ki^i or k=i. Unfortunately, this fast min-max algorithm may 
afford wrong solutions if h is strongly nonlinear on an input 
ay with an a-cut of large amplitude. However, further to the 
fact that this occurrence is supposed to be rare in practice, a 
more general algorithm would also be much more time 
consuming. On the contrary, Tan (2008) uses an even faster 
and simpler form of the above algorithm to solve the back- 
stage case, in which the signs of the partial derivatives in the 
right members of Eq. 11 are assumed a priori and so, 
Eq. 12a,b can be dispensed. The price paid is that the 
conditions of correct application of the algorithm become 
also more restrictive, as explained in Heijungs and Tan 
( 2010 ). 

Anyway, the result of applying fuzzy interval arithmetic 
to solve the uncertainty propagation problem of the back- 
stage case is the possibility distribution n for the output 
variable Y=g x shown in Fig. 5, computed for 11 a-cuts 
evenly spaced in the range of a between 0 and 1. 


7.2 Probabilistic approach for independent input variables 

A general solution for the propagation problem Y=h(X u 
X 2 , . X n ) in probability and statistics theory, even in 
the simplest case of independent input variables, is not 


available in analytical form. However, the method of 
simulation of Monte Carlo affords a numerical solution 
in three straightforward steps (e.g., Andre 2008): 

1. Generate independent TV-dimensional random samples 

A it N =(x if i, x if2 , . x iN ) for each one of the input varia¬ 
bles Xj (z=l, This step requires the availability of 

a uniform probability density generator in the interval 
[0,1], such as, the function RAND() of EXCEL™, and is 
back up on a well-known theorem of probability theory. 
More specifically, if Zi j=RAND() and X t is a continuous 
random variable with invertible distribution function F h 
then it is simply: 

Xij = Fi~ l (zi j) for i= 1,2, and j= 1,2, ...,N 

(13a) 

2. Generate an TV-dimensional sample A N =(y u y 2 , ..., y^) 
for the output variable Y, from the n samples A i>N 
0=1, n) generated in step 1, in the following 
way: 

yj = h(xij,x 2 j, ...,x nJ )fovj= 1,2, ...,N (13b) 

3. Calculate “empirically”-i.e., ultimately based on meth¬ 
od 1 of Section 4— from the sample A N generated in 
step 2, the desired probability distribution of Y (e.g., the 
probability density function f Y , if Y is a continuous 
random variable), or any of its characteristic moments 
(e.g., mean i± Y mode y 0 or standard deviation cry), with 
the help of appropriate sample estimators. 

The numerical solutions afforded in step 3 for function f Y 
or parameters (g Y yo, (J Y ) are affected by random errors with 
vanishing standard deviations for increasing sample dimen¬ 
sion N. Unfortunately, to bring the values of these standard 
deviations below acceptable thresholds, the value of N must 
usually be very large, of the order 0(TV)~ 10 4 — 10 5 , which 
turns out the method computationally heavy, especially if 
the number n of input variables is large or the functions F t rl 
(/= 1, 2, ..., n) or h are, themselves, heavy to compute. 


Fig. 5 Possibility distribution n 
(white circles denote calculation 
points) of the output variable Y 
[kg] =gi=(total amount of 
pollutant D emitted to the 
environment) that results of the 
possibilistic uncertainty analysis 
of the backstage case versus the 
possibility distributions n' (black 
circles denote calculation points) 
and 7 r" that result of the “most 
likely” probabilistic uncertainty 
analysis of the same case 
(cf. Section 7.3) 
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7.3 Comparison of possibilistic and probabilistic approaches 
for independent input variables 

The analysis of the backstage case of Section 7.1 is consid¬ 
ered illustrative of the possibilistic approach for the uncer¬ 
tainty propagation problem. For reasons that are explained 
in Section A2 of the Electronic supplementary material, the 
probabilistic uncertainty analysis that compares best with 
the former possibilistic analysis must perform the three 
following steps orderly: 

1. To convert the triangular possibility distribution tt i(x mi , 
x 0i , x Mi ) postulated for each uncertain input variable X t 
(i= 1, ..., 4) (see tt\ of Xi=a n in Fig. 6), into the 
following “most likely” probability density function 
f=U 7r f ): 



Ki{x) 


0 


da 

A X h a 


for any x E U 



In Eq. 14, \Ax h a\ designates the amplitude of the a- 
cut A Xi: a of the possibility distribution i The analytic 
solution for f is (see f\ in Fig. 6): 

1 , f ln((x 0; - x)/(x 0 . - x m .)) for x nii <x< x 0 . 

1 x Mi ~ x mi \ ln((x - x 0 ,) / (x Mi ~ x 0; )) for x 0 . < x < x Mi 

(15a, b) 

2. To apply the method of simulation of Monte Carlo 
(Section 7.2) to the propagation problem resulting from 
step 1, to obtain the probability density function f Y of the 
output variable Y. For this step to be taken, the distribu¬ 
tion function F t corresponding to f must firstly be 
obtained and, afterwards, inverted to get F t rl . The latter 



Fig. 6 Probability density function /j (Eq. 15a,b) “most likely” to the 
triangular possibility distribution 7iq of the input variable X\=a u (x m \ = 
0.5, x O i=0.6, x M \ = 1.0) of the backstage case 


has no analytic expression but can be fitted with a sixth- 
degree polynomial across the whole support \x mi , x Mi ] 
of variable X h with relative error smaller than 1% ev¬ 
erywhere except, possibly, near the extremes, where it 
may attain 3%. Eventually, the function f Y can be esti¬ 
mated with negligible random error, generating a sam¬ 
ple with dimension A=5xl0 4 and partitioning the 
relevant part of the domain of Tinto m =32 even classes. 

3. To convert the probability density function f Y obtained 
in step 2, into the “most likely” possibility distribution 
tt ,= T 2 (f Y ) (Eq. 6a, b; see Fig. 5). 

Comparing the curves of the two possibility distributions tt 
and tt' in Fig. 5, the following can be concluded: (1) the 
construction of it is much less time consuming than the one 
of tt' (in this case, tt required only 21 computations of the 
pertinent scalar function h while tt' required 5 x 10 4 ); (2) the 
modey 0 =/z(xob * 02 ? * 03 ? * 04 ) of tt is the best point estimate of 
the “correct value” y of the output variable Y (cf. Section 5.2), 
while, relative to y 0 , the mode y 0 ' of tt' is biased downwards in 
6%; and (3) tt' is much narrower than tt (the comparison is 
simplified by substituting tt' for tt" (y)=Tr' (y+%), with ^o = Yo'~ 
y (h as shown in Fig. 5), which, following the interpretation 
given in Section 5.2, means that tt' (tt") keeps much more 
uncertainty information than tt. Notice that the bias of tt' 
mentioned in conclusion 2, which is due to a bias of the 
underlying probabilistic uncertainty analysis of the type men¬ 
tioned in Section 4 (in fine), does not represent any loss of 
uncertainty information of the probabilistic approach, as the 
point-estimate y 0 of y is not an exclusive of the possibilistic 
approach. Specifically, the narrowest confidence/prediction in- 

terval /*, fory , with any desired confidence level a'<=[0,1], that 
can be obtained from (f Y y 0 ), is simply the a-cut of the trans¬ 
lated possibility distribution tt" (see above) for a=l~a'. 

Summing up, the price paid by the possibilistic approach 
for being so light-computing is the poverty of its output 
uncertainty information (the possibility distribution tt) on Y, 
which consists: (1) in the mode y 0 of tt, which is, in fact, a 
good point-estimate of the “correct value” y ; but (2) for any 
other a-cut of tt with a<\, in “confidence intervals” fory 
with confidence levels a'=l~a but amplitudes much larger 
than the optimum (compare the amplitudes of the a-cuts of 
tt and tt" in Fig. 5). On the contrary, the probability ap¬ 
proach, despite being much more time consuming, affords a 
probability distribution for Y that contains much richer out¬ 
put uncertainty information on Y. 


7.4 Probabilistic and possibilistic approaches for correlated 
input variables 

When the random input variables X t are mutually depen¬ 
dent or correlated, the step 1 of the probabilistic method 
of simulation of Monte Carlo of Section 7.2 is readily 
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adaptable to cope with it. Take, for instance, the case of 
three dependent input continuous random variables, X u 
X 2 , and X 3 . Further to the usual marginal distribution 
function F^xi) of X l9 the method now requires, as input, 
instead of the usual marginal (ID) distribution functions 
F 2 (x 2 ) of X 2 and F 3 (x 3 ) of X 3 , the conditional (respec¬ 
tively, 2D and 3D) distribution functions F 2 \ X {x 2 \x\) and 
F 3 \\^ 2 {x 3 \x { ,x 2 ), which only collapse into F 2 (x 2 ) and F 3 (x 3 ), 
respectively, when the variables X h X 2 and X 3 are mutually 
independent. Afterwards, the sample A liN is generated firstly, 
as usual, with Eq. 13a, but the samples A 2 ^ and A 3 ^ must be 
subsequently generated, the latter after the former, thus: 


X2j = F 2 \\ 1 [z 2 j\xij) and xy = F 3 \ ln2 1 (zy\x \j,x 2J ) 

(16a,b) 


So, the method of simulation of Monte Carlo has no 
limitations of principle to deal with correlated input varia¬ 
bles, and does not even become computationally heavier, 
the only problem being that it now requires more, and also 
more complex to obtain, input probabilistic information than 
in the case of independent input variables. 

On the contrary, the fuzzy interval arithmetic of Section 7.1 
not only does not require any modification to be able to cope 
with dependent possibilistic input variables but, as long as they 
depend on a common source of information (see Section A5 of 
the Electronic supplementary material and also Gebhardt and 
Kruse 1994) does not even require any more input information 
than in the case of independent input variables. 


7.5 Mixed possibilistic and probabilistic approaches 

Just for reference, Dubois (2006) mentions the works of Ferson 
and Ginzburg (1995), Guyonnet et al. (2003), and Baudrit et al. 
(2005) as recent proposals to solve the uncertainty propagation 
problem in cases in which some of the input variables are 
independent random variables and the others are possibilistic 
variables based on a common source of information, indepen¬ 
dent of the source(s) of information of the random variables. 
These authors apply a hybrid method of simulation of Monte 
Carlo (for the random variables) and fuzzy interval arithmetic 
(for the possibilistic variables) to obtain so-called upper and 
lower distribution functions for the output variable. However, 
it is acknowledged that both the (in)dependence relations 
among random and possibilistic variables, and the physical 
interpretation of the results require more research. 


8 Conclusions 

The following conclusions are highlighted: 

1. On attaching a probability measure to an uncertain 
physical variable X, care should be taken in relating 


the “correct value” x of X with parameters of the 
probability distribution of X (e.g., mean /i or mode x 0 ) 
(Section 4, in fine). 

2. Possibility and necessity measures are not indepen¬ 
dent but dual measures, i.e., they express the same 
uncertainty information in different ways (introduction 
of Section 5). 

3. In Section 5.1, three methods are proposed to convert 
available probability uncertainty information on input 
variables (e.g., more or less reliable probability distri¬ 
butions, or point and interval estimates of the “correct 
values”) into physically meaningful possibility uncer¬ 
tainty measures for the same variables, at the cost of 
losing more or less uncertainty information. On the 
construction of possibility measures directly from the 
physical non-deterministic experiments behind the un¬ 
certain estimation of input variables, there is still much 
to be done. 

4. The uncertainty information conveyed by the possibility 
distribution 7r of an output variable Y consists in (Sec¬ 
tion 5.2): (1) the mode y 0 , which is the best point 
estimate of the “correct value” y and (2) any a-cut 
(ae[0,l]), which can be roughly interpreted as a confi- 
dence interval fory with confidence level a'=\—a but 
with amplitude that is usually much larger than the 
optimum. The former interpretations are still not easy 
to extend to variables Y with multimodal possibility 
distributions n (cf. Section A2.2, Electronic supplemen¬ 
tary material), which, fortunately, are rare in practice. 

5. When it is not possible to attach a probability measure 
to an input variable X in a fully consistent way, the only 
available uncertainty information are subjective esti¬ 
mates of largely ignorant experts or it is simply wished 
to spare computation time on solving the subsequent 
uncertainty propagation problem, X can be treated as a 
possibilistic or fuzzy variable (Section 6). The second 
cause (i.e., the limitation of probability theory to deal 
with states of large ignorance) still needs to be better 
understood. 

6. Fuzzy interval arithmetic (FIA) is a computationally 
cheap way of solving the uncertainty propagation prob¬ 
lem when the input variables are fuzzy. In Section 7.1, 
very efficient but more or less general closed forms of 
FIA are given specifically to perform the uncertainty 
analysis of LCI problems. However, although not sup¬ 
posed to be restrictive in most practical cases, a deep 
physical understanding of the fundamental limitations 
of FIA and also of some more specific limitations of its 
implementations (e.g., Heijungs and Tan 2010) is still 
lacking. 

7. The classical way of solving the uncertainty propagation 
problem for random input variables through the algo¬ 
rithm of simulation of Monte Carlo can be used both for 
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independent and correlated input variables but, in the 
latter case, much more input information is required 
(Sections 7.2 and 7.4). The practical investigation of 
probabilistic correlations among input variables of 
LCI, both of their existence and of the means to get 
the additional uncertainty information required, is sup¬ 
posed to be still very incipient. 

8. On solving the uncertainty propagation problem, the 
possibilistic approach is computationally much faster 
but affords also much poorer output uncertainty infor¬ 
mation than the probabilistic approach (Section 7.3). 
However, the natural possibilistic point estimate of 
the “correct value” of an output variable [cf. Con¬ 
clusion 4, (1)] should also be adopted in the prob¬ 
abilistic approach. 

9. There exist recent proposals to solve the uncertainty 
propagation problem through mixed possibilistic/ 
probabilistic approaches when the set of input varia¬ 
bles is composed both of random and fuzzy variables, 
but the physical interpretations of their conditions of 
application and output information are still deficient 
(Section 7.5). 
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